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Perfect fluid with kinematic self-similarity is studied in 2 + 1 dimensional spacetimes with circular 
symmetry, and various exact solutions to the Einstein field equations are given. In particular, these 
include all the solutions of dust and stiff perfect fluid with self-similarity of the first kind, and all the 
solutions of perfect fluid with a linear equation of state and self-similarity of the zeroth or second 
kind. It is found that some of these solutions represent gravitational collapse, and the final state of 
£f) ■ the collapse can be either black holes or naked singularities. 
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I. INTRODUCTION 



One of the most outstanding problems in gravitation theory is the final state of a collapsing massive star, after it has 
exhausted its nuclear fuel. In spite of numerous efforts over the last three decades or so, because of the (mathematical) 
■ complexity of the problem our understanding is still limited to several conjectures, such as, the cosmic censorship 
conjecture [1] and the hoop conjecture [2]. To the former, many counter-examples have been found [3], although it is 
still not clear whether those particular examples are generic. To the latter, no counter-example has been found so far 
in four-dimensional spacetimes, but it has been shown recently that this is no longer true in five dimensions [4]. 

Lately, Brandt et al. have studied gravitational collapse of perfect fluid with kinematic self-similarities in four- 
dimensional spacetimes [5], a subject that has been recently studied intensively (for example, see [6] and references 
therein.). In this paper, we shall investigate the same problem but in 2+1 gravity [7]. The main motivation of such a 
study comes from recent investigation of critical collapse of a scalar field in 3D gravity [8-11]. It was found that in the 
3D case the corresponding problem is considerably simplied and can be studied analytically. In particular, Garfinklc 
first found a class, say, S[n], of exact solutions to Einstein-massless-scalar field equations, and then Garfinkle and 
Gundlach studied their linear perturbations and found that the solution with n = 2 has only one unstable mode [9]. 
By definition this is a critical solution, and the corresponding exponent, 7, of the black hole mass, 



M B h oc (p - p* 



is 7 = = 4/3, where k\ denotes the unstable mode. Although the exponent 7 is close to the one found 

numerically by Pretorius and Choptuik, which is 7 ~ 1.2 ±0.02 (but not the one of Husain and Olivier, 7 ~ 0.81), this 
solution is quite different from the numerical one [8]. Using different boundary conditions , Hirschmann, Wang and 
Wu found that the solution with n = 4 has only one unstable mode [11]. As first noted by Garfinkle [9], this n = 4 
solution matches extremely well with the numerical critical solution found by Pretorius and Choptuik [8]. However, 
the corresponding exponent 7 now is given by 7 = l^il -1 = 4, which is significantly different from the numerical 
ones. In this paper we do not intend to solve these problems, but study some analytical solutions that represent 
gravitational collapse of perfect fluid in 2 + 1 gravity. 
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1 Comparing the two sets of boundary conditions one will find that the only difference between them is that, in addition to 
the ones imposed by Garfinkle and Gundlach [9], Hirschmann, Wang and Wu [11] further required that no matter field should 
come out of the already formed black holes. This additional condition seems physically quite reasonable, and has been widely 
used in the studies of black hole perturbations [12]. 
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It can be shown that the conception of kinematic self- similarities given by Carter and Henriksen [13] in four- 
dimensional spacetimes can be easily generalized to D-dimensional spacctimcs with the metric 

ds 2 = I 2 { 7q6 (x c ) dx a dx b - s 2 (x c ) H i:j (x k ) dx l dx 3 } , (1.1) 

where I is an unit constant with the dimension of length, so that all the coordinates x^ and metric coefficients "f a t 
and Hij are dimcnsionlcss. Here we use lowercase Latin indices, such as, a, 6, c, ...,g, to run from to 1, lowercase 
Latin indices, such as, i, j, k, to run from 2 to D — 1, and Greek indices, such as, p, v, to run from to D — 1. 
Clearly, the above metric is invariant under the coordinate transformations 

x a = x a (x' c ) , x i = x l . (1.2) 

On the other hand, the energy-momentum tensor (EMT) for a perfect fluid is given by 

Tpv = {p + p)u^u v - pg^, (1.3) 

where w M denotes the velocity of the fluid, p and p are, respectively, its energy density and pressure. Using the 
coordinate transformations (1.2) we shall choose the coordinates such that, 

^ = (.9oo) 1/2 <5°, ,g O i-0. (1.4) 
This implies that the coordinates are comoving with the perfect fluid. Then, metric (1.1) can be cast in the form, 

ds 2 = I 2 {e 2 ^dt 2 - e 2 ^dr 2 - r 2 S 2 (t,r)H tJ (x k ) dx l dx^. (1.5) 

Following Carter and Henriksen, we define kinematic self-similarity by 

Ct V = 2^, £ 5 ti" = -era", (1.6) 

where denotes the Lie differentiation along the vector field a is a dimensionless constant, and is the project 
operator, defined by h^ v = — u^u v . When a — 0, the corresponding solutions are said to have self-similarity of 
the zeroth kind; when a = 1, they are said to have self-similarity of the first kind (or homothetic self-similarity); and 
when a ^ 0, 1, they are said to have self-similarity of the second kind. 

Applying the above definition to metric (1.5), we find that the metric coefficients must take the form, 

*(t,r) = *(*), tf(t,r) = *(aO, S(t,r) = S(x), (1.7) 
where the self-similar variable x and the vector field £ M are given by 

,u d d 9 

f M = h r — 

q dxv dt cV 

x = ln(r) -t, (a = 0), (1.8) 

for the zeroth kind, and 

d _ at d_ +r _d_ 

dx^ dt dr ' 

x = ln(r) - - ln(-i), (a ^ 0), (1.9) 
a 

for the first (a = 1) and second (a ^ 1) kind, respectively. 

As mentioned above, in this paper we shall study perfect fluid with kinematic self-similarity in 2 + 1-dimcnsional 
spacetimes with circular symmetry, for which metric (1.5) becomes 

ds 2 = I 2 {e 2 ^dt 2 - e 2 ^dr 2 - r 2 S 2 (x)d6 2 }, (1.10) 

where the hypersurfaces 9 = 0, 2tt are identified. It should be noted that the above metric is of invariance under the 
following transformations, 

t = Ai, r = Bf, g^ = C 2 g^ Vl (1.11) 
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for self-similar solutions of the first and second kinds, and 



t = i+A, r = Bf, g^ = C 2 g^, (1.12) 

for self-similar solutions of the zeroth kind, where A, B and C are arbitrary constants. 

On the other hand, to have circular symmetry, some physical and geometrical conditions needed to be imposed [14]. 
In general this is not trivial. As a matter of fact, only when the symmetry axis is free of spacetime singularity, do we 
know how to impose these conditions. Since in this paper we are mainly interested in gravitational collapse, we shall 
assume that the axis is regular at the beginning of the collapse, so that we are sure that the singularity to be formed 
later on the axis is due to the collapse of the fluid. Following [15] we impose the following conditions: 

(i) There must exist a symmetry axis, which can be expressed as 



X 



= \9ee\ -> 0, (1.13) 



as r — > + , where we have chosen the radial coordinate such that the axis is located at r — 0, and £^ is the Killing 
vector with a close orbit, given by £? e \d a = do- 

(ii) The spacetime near the symmetry axis is locally flat, which can be written as [16] 

™^-l, (1.14) 

as r — > + , where ( ) )Ce = d{ )/dx a . Note that solutions failing to satisfy this condition are sometimes acceptable. 
For example, when the left-hand side of the above equation approaches a finite constant, the singularity at r — can 
be related to a point-like particle [17]. However, since here we are mainly interested in gravitational collapse, in this 
paper we shall assume that this condition holds strictly at the beginning of the collapse. 

(hi) No closed timelike curves (CTC's). In spacetimes with circular symmetry, CTC's can be easily introduced. To 
ensure their absence we assume that 



9ee = e (e) C {6 )9^<^ (1-15) 

holds in the whole spacetime. 

In addition to these conditions, it is usually also required that the spacetime be asymptotically flat in the radial 
direction. However, since we consider solutions with self-similarity, this condition cannot be satisfied by such solutions, 
unless we restrict the validity of them only up to a maximal radius, say, r — ro(t). Then, we need to join the solutions 
with others in the region r > ro(t), which are asymptotically flat as r — ► oo. In this paper, we shall not consider such 
a possibility, and simply assume that the self-similar solutions are valid in the whole spacetime. 

It should be noted that the regularity conditions (1.13)-(1.15) are invariant under the transformations of Eq.(l.ll) 
or Eq.(1.12). Using these transformations, we shall further assume that 

$(t,0)=0, (1.16) 

that is, the timelike coordinate t measures the proper time on the axis. 

Moreover, in the analytical studies of critical collapse, one is usually first to find some particular solutions by 
imposing certain symmetries, such as, homothetic self-similarity, and then study their perturbations, because by 
definition critical solutions have one and only one unstable mode. In this paper we shall be mainly concerned with 
the first question, and leave the study of linear perturbations to another occasion. In particular, the rest of the paper 
is organized as follows: Exact solutions of the Einstein field equations with self-similarity of the first, second, and 
zeroth kinds will be given and studied, respectively, in Sees. II, III, and IV, while in Sec. V, our main conclusions 
are presented. There are also two appendices A and B, where in Appendix A the general expression of the Einstein 
tensor, among other things, for spacetimes with self-similarities of the zeroth, first and second kinds is given, while in 
Appendix B the apparent horizon is defined in terms of the self-similar variables. 



II. SELF-SIMILAR SOLUTIONS OF THE FIRST KIND 



In this section, we shall study solutions with self-similarity of the first kind. Substituting Eqs.(1.3), (1.4) and (A. 14) 
with a = 1 into the Einstein field equations G M „ = kT^ v , where k is the Einstein coupling constant, we find that 



3 



y, x - (i + y) (*, x - y) - »*, x = o, (2.1) 

'!'., .,• • $a - *,x - » - 2) 

and 

P=^-(^ +>P -* ) -^)' 

P = -^^(*^*-')-^) > (2.3) 

where y = S, x /S. Note that in writing the above expressions we had set k = 1. In the rest of this paper we shall 
continuously choose units so that this is true. Clearly, to determinate the metric coefficients completely, we need to 
impose the equation of state for the perfect fluid. In general, it takes the form [18], 

p = p(T,H), p = p(T,E), (2.4) 

where T is the temperature of the system and S the entropy. However, in some cases the system is weakly dependent 
on T, so the equation of state can be written approximately as p = p(p). In the following we shall show that in the 
latter case the only equation of state that is consistent with the symmetry of homothetic self-similarity is the one 
given by 

p=kp, (2.5) 
where k is an arbitrary constant. To show this, let us first write Eq.(2.3) in the form 



(2.6) 



P=#. (2.7) 



Then, from Eq. (2.6) we find that 

dx 

dj^p) /'(*) 

where a prime denotes the ordinary differentiation with respect to x. Inserting Eq.(2.8) into Eq.(2.7) we find that 



P-l^i. (2.9) 

which shows that in general p is a function of r and p. Taking partial derivative of the above equation with respect 
to r, and then setting it to zero, we find that 

dp{r, p) = _ 2 50f) , 1 dg(x) dx d (r 2 p) 
dr r 3 r 2 dx d(r 2 p) dr 

2/ff/V A 0) (210) 



r 3 f V 9 f 



which gives g(x) = kf(x). Then, Eq.(2.5) follows. In the following we shall further assume < k < 1 so that all the 
energy conditions hold [19] and that the pressure is non-negative. 
The combination of Eqs.(2.3) and (2.5) immediately yields 



»— rb- (2 - n) 



Then, Eqs.(2.1) and (2.11) have the solutions, 



S(x) =S e- x / (1+k \ (2.12) 
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while Eq.(2.2) becomes 



+ k (*, x + -j^jj [(* - 1) *,* - 1] 



_ e 2 [xfT + (l-fe)*-*o] 



where $o and So are integration constants. 



* xx + (1 - fc) * X ( * X + 



fc + 1 



= 0, 



(2.13) 



A. Stiff Fluid (fc = 1) 

When fc = 1, i.e., the stiff fluid, Eq.(2.13) has the general solution, 

*{x) = gin (1 - e*"* °) - ± (a: - x ) + *o, 
$(x) = g ln(l-e x - x °) + x + * , 
S(x) - Soe^ 2 , 



(2.14) 



where ^0: £o> So and 9 are the integration constants, and xq = 2$o. Imposing the regularity conditions of Eqs.(1.13)- 
(1.15) and the gauge condition (1.16), we find that 



So = 2e-*°/ 2 ) *o 



-xo, 



(2.15) 



while xo and q are arbitrary and cannot be removed by the transformations of Eq.(l.ll). Then, the general solution 
is given by 



ds 1 



{0 



I 2 1 



3 ) 



2? 



(x+x ) rfr 2 _4 r 2 e -(x+xo) rf6/ 2| 



and the corresponding energy density of the fluid reads, 



P = P 



1 - 2q 



4Z 2 (-t) 2 (1 - e x ~ x °) 



2<r 



(2.16) 



(2.17) 



It should be noted that in writing Eq.(2.16), we had implicitly assumed that x — x < 0. On the other hand, from 
Eq.(2.17) we can see that we must have q < 1/2 in order to have p > 0, a condition that we shall assume in the rest 
of the paper. The metric is singular on the hypersurface x — x , and depending on the values of q, the nature of the 
singularity is different. 



1. < q < 1/2 

In this case the singularity is a curvature one and marginally naked. In fact, from Eq.(B.lO) we find that 



7; = 



fe x 



2l 2 r (1 - e x - x °) q 
ge x 



,(x -x)/2 



- 1 



2Z 2 r(l 



,(x -x)/2 _|_ j 



(2.18) 



from which we can see that 6 n < for any given x G [0,xo], and 8i is positive for x < xo, zero only on the surface 
x = x . Thus, now the hypersurface x = x is a marginally trapped surface and the spacetime singularity located 
there is marginally naked. On the other hand, for any given hypersurface x = Const., say, C, its normal vector is 
given by 



(2.19) 
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from which wc find 



n a npg alj = - 




(1 




l-2q 



(2.20) 



Clearly, on the hypersurface x — x the normal vector, n a , becomes null. That is, the nature of the singularity is 
null. The corresponding Penrose diagram is given by Fig. 1. 



FIG. 1. The Penrose diagram for the solutions given by Eq.(2.16) with < q < 1/2. The spacetime is singular on the double 
line x — xo, which is a null surface and on which the expansion of the outgoing radial null geodesies is zero, #;(£, r)\ x=xa = 0, 
but 9 n (t,r) < in the whole spacetime, including the hypersurface x — xo- 



In this case the metric is free of spacetime and coordinate singularities on the hypersurface x — x , and valid in the 
whole region t < 0, r > 0. However, now the spacetime becomes singular on the hypersurface t — [cf. Eq.(2.17)]. 
From Eq.(2.18), which is also valid for q = 0, we can see that this singularity is not naked, and always covered by the 
apparent horizon formed on the hypersurface x — xq. In the region x > xq, which is denoted as Region / in Fig. 2, 
Eq.(2.18) shows that 9i becomes negative and 6i6 n > 0, that is, all the rings of constant t and r now become trapped. 
Thus, in this case Region / can be considered as the interior of a black hole, which is formed by the gravitational 
collapse of the stiff fluid in Region //. The corresponding Penrose diagram is given by Fig. 2. 



o 




2. q = 



t = 



FIG. 2. The Penrose diagram for the solutions given by Eq.(2.16) with q = 0. The spacetime is singular on the double line 
t = 0. All the rings of constant t and r are trapped in Region 7 but not in Region II, where I = {x a : x > xo,t < 0, r > 0} 
and II = {x a : x < xo, t < 0, r > 0}. The hypersurface x = xo is a null surface and represents an apparent horizon. 



It should be noted that apparent horizons and black holes are usually defined in asymptotically flat spacetimes [19]. 
As mentioned previously, the spacetimes considered here are not asymptotically flat, therefore, strictly speaking the 
above definitions are, in some sense, the generalization of those given in [19]. As a matter of fact, to be distinguishable 
with the asymptotical case, Hay ward called such apparent horizons as trapping horizons and defined black holes by 
future outer trapping horizons [20]. For the sake of simplicity and without causing any confusions, following [9-11] 
we shall continuously use the notions, apparent horizons. 



3. q < 

In this case the spacetime is free of curvature singularity on the hypersurface x — x , although we do have a 
coordinate one there. Then, to have a geodesically maximal spacetime, we need to extend it beyond this surface. 
Introducing two null coordinates u and v via the relations 

t = ~\ [(-«)" + (-v) n f, 

r=l e *° [(-«)" -(-t,)"] 2 , (2.21) 



where 



(2.22) 



2^ + 1' 

we find that the metric (2.16) becomes 

ds 2 = I 2 ^2e 2a ^ v) 'dudv - R 2 (u,v)d6 2 } , (2.23) 

with 

a(u, v) = (1 - 2q) In [(-«)" + (-«)"] + a , 
R(u,v) = {-u) 2n - (-v) 2n , 

<7 = \ In (2n 2 4 2 «) , (2.24) 
and the corresponding velocity and energy density of the stiff fluid are given, respectively, by 

Po (uv) n - 1 _ _ 2n(2n-l)^_ 2ao 



Po = ^l -e- 2ao - (2.25) 



From Eq.(2.21) we can see that the region t < 0, r > 0, x < xq in the (t, r)-plane has been mapped into the region 
u < 0, v < and v > u in the (u,w)-plane, which will be referred as to region II, as shown in Fig. 3. The null 
hypersurface x = x , as can be seen from Eq.(2.20), is mapped to the one v = 0. Region I, where u < 0, v > 0, |u| > v, 
represents an extended region. In this extended region, the metric is real only for n ^ (2m + l)/(2j), where j and 
m are integers. When n = (2m + l)/(2j), a possible extension may be given by replacing —v by \v\. However, this 
extension is not analytical. In fact, no analytical extension exists in this case. The only case where the extension is 
analytical is the one where n is an integer. When the extension is not analytical, it is also not unique. Thus, to have 
an unique extension in the following we shall assume that n is an integer. On the other hand, from Eq.(B.7) we find 
that 

Oi = 2j ^ L (-v? n -\ On - - 2j ^(-uf n -\ (2-26) 
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where R = rS. From the above expressions we can see that 9 n is always negative in both regions I and II, while 9i is 
positive in Region II, zero on the hypersurface v = 0, and negative in the extended region, I. Near the surface v = 
the only non-vanishing component of the energy-momentum tensor is given by 



T — 

-L 11.11. 



l 2 p e 



2a 



{-Uf 



(2.27) 



which represents the energy flow moving from Region II to Region / along the hypersurfaces u — Const. To study 
the above solutions further, let us consider the two cases n = 2m + 1 and n = 2m, separately, where m = 1, 2, 3, .... 

When n is an odd integer, i.e., n = 2m + 1, from Eq.(2.25) we can see that the spacetime is singular on the 
hypersurface u = —v in the extended region /, on which we have R(u, v) = 0. On the other hand, all the rings of 
constant t and r are trapped in this region, as we can see from Eq.(2.26), now we have 9i(v > 0) < and 9i9 n > 0. 
Thus, region I can be considered as the interior of the black hole that is formed from the gravitational collapse of the 
fluid in region II. The corresponding Penrose diagram is that of Fig. 3. 



R = 




FIG. 3. The Penrose diagram for the solutions given by Eqs. (2.23) and (2.24) with n being an integer. The rings of constant 
t and r are closed trapped surfaces in Region 7, but not in region II. The dashed line v = 0, or x — xo, represents an apparent 
horizon. When n is an odd integer, the spacetime is singular on the horizontal double line R = 0, while when n is an even 
integer, the spacetime has an angular defect there. 



When n is an even integer, say, n — 2m, from Eq.(2.25) we can see that no spacetime curvature singularity is 
developed on the axis R = in the extended region I, although all the rings of constant t and r are also trapped in 
this region, as now we still have $i(v > 0) < and 9i(v > 0)9 n (v > 0) > 0. However, the local flatness condition 
(1.14) is not satisfied there. In fact, it can be shown that now we have 



AX 



a/3 



+ 1, 



(2.28) 



as v — > — u. Thus, unlike that on the axis R = in Region II, where the local flatness condition is satisfied, now 
the spacetime on the axis R = in region / has angle defect. The corresponding Penrose diagram is also given by 
Fig. 3, but now the horizontal double line R = 0, instead of representing a curvature singularity, now represents an 
angle-defect-like singularity [17]. 

Moreover, in the extended region I, where uv < 0, the function p becomes negative, and the three-velocity 
becomes imaginary, as we can see from Eq.(2.25). A close investigation shows that the energy-momentum tensor in 
this extended region actually takes the form 



= p(2r tt r v + g^) , 



(2.29) 



where 



le a ( 1 

^ = T2 _ ~ s; 



Pa \uv\ 



2m-l 



(2m-l)/2 



( u 2m _|_ v 2m)( 6m - 1 )/ m 



, (n = 2m), 



(2.30) 
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from which wc find that 



gTr^ = -1. (2.31) 
Thus, in the extended region the source is no more a stiff fluid. Introducing the unit vectors, 

^ = ^^" + 1^, 9, = lR5l (2.32) 

where u v u v = 1 and Q V Q V = —1, we find that Eq.(2.29) reads 

Tft V = pu^u v - p r r ll r v - Pe6^6 v , 

P = -Pr =Po = P ° lUVl2 V 1)f ■ (2.33) 

Thus, the source in region I now becomes an anisotropic fluid with its energy density p and two principal pressures, 
p r and pg, in the direction of , respectively, and 9„. Note that although the pressure in the r„ direction is negative, 
the anisotropic fluid satisfies all the three energy conditions [19]. With this odd feature, it is not clear whether the 
corresponding solution can be interpreted as representing gravitational collapse of the stiff fluid. In fact, if we consider 
the fluid from its microscopic point of view, we might be able to rule out such a change in the equation of state across 
the hypcrsurfacc v = [21], a problem that is under our current investigation. 

Before turning to study other cases, we would like to note that it is well-known that a still perfect fluid is energetically 
equivalent to a massless scalar field when the velocity of the fluid is irrotational V^u^ = and the massless scalar 
field is timelike (p ta (f)^g a ^ > 0. Clearly, these conditions are fulfilled in Region 77. In fact, comparing the solutions 
given by Eq.(2.24) with the corresponding massless scalar field ones found in [11], one can find that they are actually 
the same in Region II. However, as shown lately in 4-dimcnsional case [22], the spacetime across the horizon v = 
can be quite different. Our analysis given above and the one given in [9,11] show that this is also the case in 2 + 1 
Gravity. In [22], it was also shown that the linear perturbations in these two cases are different. Thus, it would be 
very interesting to study the linear perturbations of the above solutions in terms of stiff fluid. 



B. Dust Fluid (k = 0) 

When k = 0, the pressure of the perfect fluid vanishes, i.e., the dust fluid. Then, it can be shown that in this case 
the general solution is given by 

<$>(x) = $ , 

tf(ar) = In (1 - e*-* ) -x + * , 

S(x) = S a e- X , (2.34) 

for which we find 

p = n-ty(i-e*-*°y p = °- (2 - 35) 

Since now we have R(t,r) — rS(x) — —Sot, we can see that these solutions are Kantowski- Sachs like [16], and may 
be interpreted as representing cosmological models. Because in this paper we are mainly interested in gravitational 
collapse, in the following we shall not consider this case in any more details. 

C. Perfect Fluid With k ^ 0,1 

When k ^ 0, 1, introducing the function Z{x) by 



Z{x) = cxp|2 
we find that Eq.(2.13) can be written in the form 



1 + k 



(1 - k) * - $ 



(2.36) 
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2Z (fc - Z) Z tXX + {Z- 3k) Z, 2 + 2Z 2 Z x 



^Z 2 (fc 2 - Z) = 0. (2.37) 



(fc+1) 

We are not able to find the general solution of this equation, but a particular one given by 

®(x) = $ , 

T , ^ a; $o 

*0) = "I 7 + 



fc+1 1-fc' 

S(x) = S e-T^r. (2.38) 



Then, it can be shown that the conditions of Eqs.(1.13)-(1.16) are fulfilled, provided that 



$o = 0, 5 , = -i-, fc>0. (2.39) 



Thus, the corresponding metric takes the form 

ds 2 = I 2 \dt 2 - e- 2 ^ 1+fe ) 



(2.40) 



and the pressure and energy density of the fluid are given by 



fc 

Z 2 (l + fc)" 2 (-t) 2 



P = k P = ;2/1 , u , 2 , 2 . (2.41) 



From the above expression we can see that the spacetime is always singular at t — 0. However, this singularity is not 
naked. In fact, from Eq.(B.lO) we find that 



_ /e^+fc) / kx/(1+k) \ 



x/(l+k) 
a _ ^ c 

P(l + fc)r 

from which we can see that n is always negative, and 



/ 2 (l + fc)r 

(k + e kx/(i+k)y ( 2 42 ) 



where xq is given by 



{> 0, x < Xq, 
= 0, x^xo, (2.43) 
< 0, X > Xq, 



Xo =l±^\n(k). (2.44) 



Thus, the singularity at t = is always covered by the apparent horizon localized on the hypersurface x — xq. In the 
region where t < and x > xo, the rings of constant t and r are all trapped, as now we have 9i < and 9i9 n > 0. 
Therefore, in the present case the collapse of the fluid always forms black holes. Moreover, the normal vector n a to 
the hypersurface x = x is still given by Eq.(2.19) but now with 

2x /{l+k) 

n a n p g af} = (l - fc 2 ) < 0, (2.45) 

since now we have < fc < 1. That is, the apparent horizon is always timelike and the corresponding Penrose diagram 
is given by Fig. 4. 
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FIG. 4. The Penrose diagram for the solutions given by Eq.(2.40) for < k < 1. The apparent horizon x = xo, denoted by 
the curve AH is always timelike. In Region II we have 9i > and 9i6 n < 0, but in Region I all the rings of constant r and t 
are trapped as now we have 6i < and 9i6 n > 0. The spacetime is singular on the double line t — 0. 



III. SELF-SIMILAR SOLUTIONS OF THE SECOND KIND 



When a ^ 1, the term that is proportional to r~ 2 has different power-dependence on r from the term that is 
proportional to i~ 2 , when they are written out in terms of r and x, since t — — r a e~ ax . Then, it can be shown that 
the Einstein field equations in this case become 

y, x - (1 + y) (*, x - y) - y$, x = 0, (3.1) 
$ ** + $ * - *,* - y - 2) = 0, (3.2) 
f,xx - * * - * * + y + 1 - a) + (1 - a) y = 0, (3.3) 



and 

















[ r 2 



" W e_2$ [(1 + V) _ (1 _ a) 2/1 } ' (3 ' 4) 

where in writing Eqs.(3.2) - (3.4), we had used Eq.(3.1). From the above equations we can see that now the Einstein 
field equations are already sufficient to determinate completely the metric coefficients ^(x) and S(x). Following 

Maeda et al [24] , it can be shown that the symmetry of the self-similarity of the second kind is inconsistent with a 
polytropic equation of the kind, 

p = kpP, (3.5) 



or 

p = kn f3 1 p = m b n + V _ ^ , (3.6) 

unless k = 1, where the constant nib denotes the mean baryon mass, and n(t,r) the baryon number density [18]. 
Thus, in the following we shall consider only the case (3—1. Combining it with Eq.(3.4), we find that 

[l + (l + k)y]*, x = 0, (3.7) 
[l + (l + fc)y]* ia! -(l-a)y = 0. (3.8) 

Since a ^ 1, we must have $ x = 0. Therefore, for the perfect fluid with the equation of state p = kp and self- similarity 
of the second kind, it must move along the radial timelike geodesies. To solve Eqs.(3.1)-(3.3) and (3.8), let us consider 
the two cases y ^ — 1 and y = — 1 separately. 
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Case A) y ^ — 1: In this case, from Eq.(3.1) we find that 



y,x 
i + y 



while from Eq.(3.8) we obtain 



[(1 - a) - (1 + fc) tt x ] . 



l + (l + fc)y 

Inserting the above expressions into Eq.(3.3) we find that it has three different solutions, 

(a) fc = 0, 
(6) = 0, 

(c) i/,x - (a - 2)y(l + 1/) = 0. 
When fc = 0, it can be shown that Eqs.(3.1)-(3.3), (3.7) and (3.8) have the general solution, 



(3-9) 



(3.10) 



(3.11) 



$(a0 = *o, *(x)=In 1 + (a - l)e a(x °- x) 



+ *o, 



1 



o a(x -x) 



(fc = 0). 



S(z) = 5 

When y,rr = 0, we find that the general solution is given by 

$(z) = $o, S(x)=S e ax , 
if!(x) = ax + <J>o, a 



(3.12) 



1 + fc' 



(3.13) 



When 7/ j:E — (a — 2)y(l + y) = 0, it can be shown that the corresponding solution is given by Eq.(3.13) with a = 2. 
Thus, this case is a particular one of Eq.(3.13). 

Case B) y = —1: In this case from Eq.(3.8) we find that 



l-a 



k 



(3.14) 



Inserting it into Eq.(3.3) we find that there are only two solutions given, respectively, by 



(i) fc = 1, 

(ii) a = 1 + fc. 



When fc = 1, the general solution is given by 



and when a = 1 + fc, we have 



$(x) = $ , 5(x) = S* e-, 
tf(a;) = (l-a)a; + tfo, ( fc = 1). 



$(a0 = * o , S(z) = S e-*, 
*(or) = -x + * , a = 1 + fc. 



(3.15) 



(3.16) 



(3.17) 



Thus, the most general solutions with self-similarity of the second kind (a ^ 0, 1) and the equation of state p = kp 
consist of four classes of solutions, given, respectively, by Eqs.(3.12), (3.13), (3.16) and (3.17). In the following, let us 
consider them separately. 
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A. y^-1, fc = 



In this case, the solutions are given by Eq.(3.12). Applying the conditions (1.13)-(1.16) to them, we find that $ = 0, 
So = e*° and a < 1. Then, using the transformations (1.11) we can further set = 0. Thus, the corresponding 
metric finally takes the form, 



ds- 



= l 2 {dt 2 ~ 



1 - (1 - a)e a ^°-^ 



n2 



dr z - r 2 



1 — e"' 10-1 ' 



d0 2 \, (a< 1), 



and the energy density and pressure of the fluid are given by 



! (1 - a)e 2a( - x °-'^ 



l 2 {at) 2 |l - e a(*o-*)| 1 1 — (1 — a)e a( - x »- x ^\ ' 
p = 0, (a < 1). 

From the above expressions we can see that the spacetime is singular on the hypersurfaces, 

(i) t — 0, (ii) x — x , (Hi) x = x\, 

where 



x\ = x In 

a 



1 



1 - a 



< xo- 



On the other hand, from Eq.(B.lO) we find that 

/ 



p r | ^ _ e u(x -x) I 

9 

Pr 1 1 — e a i x a—x) I 



1+ — 



(0 < a < 1), 



for < a < 1, and 



01 = 



f 



p r | ^ _ e a(x -x) | 

9 

\2 r I ^ _ e a{xa—x) I 



1 + 



1- ( — 



(a < 0), 



for a < 0, where ta = e 



-ai /(l-a) 



(3.18) 



(3.19) 
(3.20) 

(3.21) 



(3.22) 



(3.23) 



From Eq.(3.22) we can see that, in the case < a < 1, 6 n is always negative for any r G [0, oo), but 0; is negative 
when r > ta, zero when r = ta and positive when r < ta [cf. Fig. 5]. That is, the spacetime is closed far away from 
the axis even at the very beginning (t = — oo). This property makes the solution very difficult to be considered as 
representing gravitational collapse. 

When a < 0, from Eq.(3.23) we can see that 9i now is always positive, but 8 n changes signs at r = ta- In particular, 
it shows that the ingoing radial null geodesies become expanding in the region r > ta- With this odd property, it is 
also very difficult to consider the corresponding solutions as representing gravitational collapse. 
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FIG. 5. The spacetime in the (t, r)-plane for the solutions given by the metric (3.18). It is singular on the hypersurfaces 
t = 0, x — xo and x = xi, where xq > x\. 



B. y = Const.y= -1) 

In this case the solutions are given by Eq.(3.13). It can be shown that the regularity conditions (1.13)-(1.15) and 
the gauge one (1.16) require <J> = 0, So = e*°/(l + a) and a < 1 + k. On the other hand, using the transformations 
(1.11), we can further set ^ = 0. Then, we find that 



ds 2 = I 2 |di 2 - e 2ax 
and the corresponding physical quantities are 

p = kp 



dr 2 + 



ka 2 



1 + a 



d6' 



, (a<l + k), 



(3.24) 



Z 2 (l + k) 2 {at) 2 



fe~ 



P(l + k)r 

g e -ax 

~l 2 (l + k)r 



(l + *0(l + o)- (Z^WTfe) 



(l + A)(l + o) + 



U+a 



(_t)fe/(l + fe) 



(3.25) 



From the above we can see that the spacetime is singular on the hypcrsurfacc t — 0, and the singularity is not naked. 
In fact, it is covered by the apparent horizon located on the hypersurface, 

t ah (t) = [(1 + + a)(-i) fe/(1+fe) ] V(1+a) . (3.26) 
The normal vector to the hypersurface r — rAii(t) = is given by n a = d(r — rAH{t))/dx a . Then, we find that 

n a n a = -I' 2 (1 + k - a ) 2 «/( 1 + fe -«) (_ t )-2(l-a)/(l+fc-a) ^ _ fc 2) _ ( g 2?) 

Thus, the apparent horizon is timelike for < k < 1, and the corresponding Penrose diagram is that of Fig. 4. When 
k = 1, it is null, and the corresponding Penrose diagram is that of Fig. 2. 



C. y = -1, k = l 

In this case the solutions are given by Eq.(3.16). The transformations (1.11) and the gauge condition (1.16) enable 
us to set $o = = 0, Sq = 1. Then, the metric takes the form 
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ds' 



= I 2 [dt 2 - di* - {-t) 2 ' a d6 2 } , (3.28) 



from which we can see that the solutions are Kantowski-Sachs like and do not satisfy the regular conditions (1.13) 
and (1.14). The corresponding energy density is given by 

>=*=iw- (3 - 29) 

Thus, to have the energy density be positive, we must assume that a > 1. 

D. y = -1, a = 1 + k 

In this case, the solutions are given by Eq.(3.17). Using the transformations (1.11) and the gauge condition (1.16) 
we can set $o = an d Sq = 1. Then, the metric takes the form, 



ds 2 



I 2 {dt 2 - e 2{x °- x) dr 2 - {-t) 2 ' a d6 2 } , (3.30) 



where a = 1 + k and x = *o- Thus, in this case the solutions are also Kantowski-Sachs like, and the corresponding 
energy density and pressure of the fluid are given by 

p = kp= ¥W- (3 - 31) 

As we mentioned in the last section, when k = 1 the perfect fluid is energetically equivalent to a massless scalar 
field. Indeed, the solutions given by metric (3.28) and the ones given by Eqs.(3.24) and (3.30) with k = 1 are the 
solutions found lately for a massless scalar field in [23] . 

IV. SELF-SIMILAR SOLUTIONS OF THE ZEROTH KIND 

In this case, combining Eqs.(1.3) and (A. 8), we find that the Einstein field equations can be cast in the forms, 

y,x - (1 + y) (*,*- V) -»*,* = 0, (4.1) 
+ - *, x - y - 2) = 0, (4.2) 

- *,x (*,x - *,x + y + 1) + y = o, (4.3) 

and 

P=^{r a e-»^-c-»*„}. 

P = ^2 {(! + y)e" 2 **,x - r 2 ^ 2 * [(1 + 2/) *,x - y]} , (4.4) 

where in writing Eqs.(4.2) - (4.4), we had used Eq.(4.1). Similar to the solutions with self-similarity of the second 
kind, now the Einstein field equations already determinate completely the metric coefficients. It can be also shown 
that the self-similarity of the zeroth kind is inconsistent with the equation of state given by Eq.(3.5), unless (3=1. 
Then, following the last section, let us consider solutions that satisfy the equation of state (2.5), which, together with 
Eq.(4.4), yields, 

[l + (l + fc)y]$, x =0, (4.5) 
[l + (l + fc)y]*, x -y = 0. (4.6) 



The above equations have the solution 

•(■>-•»• *-* = iT(rW (4 - 7) 

Thus, similar to the case of self-similarity of the second kind, now the fluid must also move along radial timelike 
geodesies. Inserting the above expressions into Eqs.(4.1)-(4.3), we find that solutions exist only for three special cases, 
k = 0, ±1. 
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A. k = 



When k = 0, the general solution is given by, 



$(x) = $ , S(x) = So (x - x) , 
tf(ar) = In (ar - a: - 1) + tf , (fc = 0). 



(4- 



Then, it can be shown that the conditions (1.13)-(1.16) require $0 = and So = e*°. On the other hand, using the 
transformations (1.12), we can further set ^o = 0. Then, the solutions are finally given by 



ds 2 = I 2 \dt 2 - (x - x - if dr 2 - r 2 (x Q - x) 2 d6 2 } , 
and the corresponding energy density is given by 



1 



I 2 (xq — x — 1) (xq — x) 



, p = 0. 



(4.9) 



(4.10) 



From this expression we can see that the spacetime is singular on the hypersurfaces x = x — 1 and x = x [cf. Fig. 
6]. From Eq.(B.9) we also find that 



0i 



f 



1 2 (xq — x)r 
9 



1 2 (xq — x)r 



(1 + r), 
(1-r). 



(4.11) 



Thus, 9i is positive for x < xq and negative for x > xo, while the signs of 9 n get changed when across the hypersurfaces 
r = 1 and x = xq. With this property, one can see that the solutions cannot be considered as representing gravitational 
collapse. 

















/ " / 






\ / 









1 




\ 1 *° 






/ / " 






\ 1* 





FIG. 6. The spacetime in the (t,r)-plane for the solutions given by Eq.(4.9). It is singular on the hypersurfaces x = xo and 
x = xo — 1. The signs of 9 n change when across the hypersurfaces r = 1 and x — xo- 



B. k = -1 

When k = — 1, the general solution is given by 

$(x) = $o, *(x)=aa; + *o, - S e ax , (4.12) 

for which we have 

P = -P=^^. (4.13) 
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This is a 3-dimensional de-Sitter solution. In fact, setting 



t = e- <Po i+— , r=[(l + o) 



r 



a 

0*0 



/^ae"* , (4.14) 



(1 + a)S 

we find that the corresponding metric can be cast in the form 

ds 2 = I 2 {dP - e- 2 ^ (df 2 + r 2 d9 2 ) } . (4.15) 



C. Jfe = 1 

In this case, the general solution is given by 

*(a0 = *o, *{x)=x + * , S{x)=S e- x , (4.16) 

for which we have 

e -2*o 

P = P = —■ (4-17) 

Clearly, now none of the three energy conditions is satisfied, and the physics of the solution is unclear (if there is any). 

V. CONCLUSIONS 

In this paper we have first generalized the notion of kinematic self-similarity of four-dimensional spacetimes to any 
dimensions for the metric given by Eq.(l.l), and then restricted ourselves to (2 + l)-dimensional spacetimes with 
circular symmetry. 

In Sec. II, we have studied solutions of the Einstein field equations with homothetic (the first kind) self-similarity 
for a perfect fluid. It has been shown that the only equation of state that takes the form p = p(p) and is consistent with 
the self-similarity is p = kp, where k is a constant. In the latter case, a master equation has been found, Eq.(2.13). 
Then, the general solutions for the stiff fluid (fc = 1) have been given in closed form. It has been found that some of 
these solutions represent formation of marginally naked singularities and the others represent the formation of black 
holes from gravitational collapse of the self-similar fluid. All the solutions of dust fluid (k = 0) have also been given, 
and found that they are all Kantowski-Sachas-like. When k 7^ 0,1, only particular solutions have been found and 
shown that the corresponding collapse always forms black holes. 

In Sees. Ill and IV, all the solutions with self-similarity of the second or zeroth kind and the equation of state 
p = kp have been found and studied. In particular, it has been shown that the fluid must move along timclikc radial 
geodesies. It has been also shown that some of those solutions represent gravitational collapse of the perfect fluid, 
while the others don't. In the collapsing case, black holes are always formed. 

In critical Type II collapse, all the critical solutions found so far have either discrete self-similarity or kinematic 
self-similarity of the first kind, and no critical solutions with kinematic self-similarity of the zeroth or second kind 
have been found, yet [25,26]. However, in Statistical Mechanics, critical solutions with kinematic self-similarity of the 
second kind seem more generic than those of the first kind [27]. Thus, it would be very interesting to study the linear 
perturbations of the solutions found in this paper with self-similarity of the zeroth or second kind. 

Finally, we would like to note that Ida recently showed that a (2+l)-dimcnsional gravity theory which satisfies the 
dominant energy condition forbids the existence of black holes [28]. This result does not contradict with the ones 
obtained here, as in our case the surfaces of apparent horizons of the black holes are not compact, a precondition that 
was assumed in [28] . 



17 



APPEND A: CIRCULARLY SYMMETRIC SPACETIMES WITH KINEMATIC SELF-SIMILARITY 



The general metric of (2 + l)-dimensional spacetimes with circular symmetry can be cast in the form, 

ds 2 = I 2 [e 2 ^dt 2 - e™^ r Ur 2 - r 2 S 2 (t, r) d6 2 } , (A.l) 

where I is a constant and has dimension of length. Then, it is easy to show that the coordinates {x^} — {t, r, 9}, 
the Christoffcl symbols, r*„, the Riemann tensor, R°\, the Ricci tensor, R^, and the Einstein tensor, G^, are all 
dimensionless, while the Ricci scalar, R 1 has the dimension of l~ 2 , and the Kretschmann scalar, TZ = R alx R^^x, 
has the dimension of Z~ 4 . 

For the metric (A.l), we find that the non-vanishing Christoffel symbols are given by 

r° = $ t , r° 01 = $ r , r° n = e 2 ^-*H, t , r° 22 = r 2 e - 2 *ss, t , 

Tl a = e 2 (*-*)$, r , rj 1= * t , = * jrs rl 2 =re-™S(S + rS, r ), 

r 2 - ^ r 2 - - + rS ' r f a 9) 

2 ~ S ' ~~ rS ' 

and the Einstein tensor has the following non-zero components 

-2* 

G tt = — {re 2 ^S. t ^. t - e 2 * [rS. rr + 2S, r - (S + rS, r ) r ]} , 

G tr = Ktr - S ,* r - !) - ( r5 ,r + S ) *■*] . 

TO 

-2* 

G 09 = -r 2 S 2 {e- 2 * [*,„ + (*, t - *, t ) *, t ] 

- e- 2 *[*, rr + ($ r -*, r )$ r ]}. (A.3) 



A.l Spacetimes with Self-Similarity of the Zeroth Kind 

To study solutions with self-similarity of the zeroth kind, let us first introduce the self-similar variables, x and r 
via the relations, 

x = ln(r) -t, T = t, (A.4) 

or inversely 

t = T, r = e x+T . (A.5) 

Then, for any given function f(t, r) we have 

J,t = f,T ~ J,xi f,r = — J.xi 

r 

f,tr — ~ (f.xx f,Tx) j f.rr = (/:XX f,x) j 

/,« - /,TT - 2/, TI + /,XX- (A.6) 

Substituting these equations into Eq.(A.3), we find that 

-2* 

G tt = j- {e 2 * [S xx + S, x - tf x (5, x + 5)] 

-r 2 e 2 * (5, T - S x ) (*, T - *, x )} , 
G tr = [S xx - S, TX + (S, x + S) (*, T - * x ) 
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+ (S, T - S, x ) - 1)] , 

-2* 

°rr = { e ™*' X {S ' X + S) ~ r2g2 * [S ' TT ~ 2S ' TX + S ' XX 

+ (S >T -S >X )(* X -* T )]}, 

Gee = S 2 e~ 2 ^ {e 2 * [$ xx + ($ x - *, x - 1)] 
-r 2 e 2 * [*, xx + *, TT - 2tf TX 

- (*, T - *,*) (*,t - *,t - $ x + *,*)]} • (A.7) 

For the self-similar solutions, the metric coefficients \P and S are functions of x only. Then, Eq.(A.7) reduce to, 

-2* 

G« = -— {e 2 * [</,x + (V + 1) (» - *,*)] - r 2 e™V, x y} , 

Gtr = J [y,x + (y + 1) {y - * ,*) - 1/*,*] , 

-2* 

G r r = {e 2 **,, (2/ + 1) - r 2 e 2 * [y, x + y (y - $ x )]} , 

G ee - ^ e -2(* + *) | g2 * [$ ^ + $ x { $ x _ y x _ 1}] 

-rV* [* xx - * x ($ x - *, x )] } , (A.8) 

where 



A. 2 Spacetimes with Self-Similarity of the First and Second Kinds 

To study these kinds of self-similar solutions, let us introduce the self-similar variables, x and t by 

x = In ( T ; I , r=-ln(-t), (A.10) 

V (-*)-/ 

or inversely, 

r = e (as-r)/a t = _ e -^ (A.H) 

where a is a dimensionless constant. When a = 1, the corresponding spacetimes are said to have self-similarity of 
the first kind or homothetic self-similarity. Otherwise, they are said to have self-similarity of the second kind. 
For any given function / (i, r), now we have 

f.t = — \ {af, T + f.x) , f, r = -f, x , 
at r 

f.tr ; \(%f.TX ~t~ f,xx) 7 f.rr — K {f.xx f.x) j 

air r z 

= ^2 ("V.rr + 2c */>™ + + ^f.r + «/,*) ■ (A. 12) 

Substituting these equations into Eq.(A.3), we find that 



Se 2 * 



r 2 

2 

_^_ e 2 * (a*, T S, T + 9 >X S, T + *, T S, X ) [> , 
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Gtr = — ^ {S, xx - *, x (S. x + S)- S, x - 1) 
+a [S, rx - * T (S, + 5) - 5 r (<& x - 1)]} , 

Gn- = ^ ? ^ 5 {a a e»[*„(S 1 , + S)] 

~~£2 e {S,xx ~ S yX <&, x ) 
2 

-^e 2 * {aS ;TT + 2S, rx - S, x (*, T - 1) 
-5, T [a (*, T -1) + *,*]}}, 
G ee = {a 2 e" 2 * [$ xx + <& x (<& x - *, x - 1)] 

-^e' 2 * [* xx - x ($ x - *, x - a)] 

2 

_^ e -2* {a y TT + 2 * rx _ * r [ a ($ r _ $ r _ i) + $ x _ $ iX ] 

-*,x(*,r-*,r)}}. (A.13) 

For the self-similar solutions, the metric coefficients are also functions of x only, but now with x being given by 
Eq.(1.9). Thus, for the self-similar solutions, Eq.(A.13) reduces to 

G « = ^2*.*» - ^e 2 <*-*> + (v + !) (» - *.*)] - 

GW = 4- [y, x + (y + 1) {y - *, x ) - y$ x ] , 

<?rr = -^e 2 <*-*) [y, x + y (y - $ x + a)] + 1$ ;X (y + 1) , 

Gee = S 2 {e- 2 * [<&** + - *,* - 1)] 



r 



2 



a 2 t 2 



e 



-2$ 



[*, xx - *, x (*, x - *, x - «)] \, (A.14) 



where y is given by Eq.(A.9). 



APPEND B: APPARENT HORIZONS OF CIRCULARLY SYMMETRIC SPACETIMES 

In [11], the ingoing and outgoing radial null geodesies were studied in double null coordinates, and apparent horizon 
was defined as the outmost hypersurface where the expansion of outgoing null geodesies vanishes. In this appendix, 
we shall write down the expansions in terms of the self-similar variables. To do this, let us first introduce two null 
coordinates u and v via the relations 

du = f (e*di - e^dr) , dv = g (e*dfc + e^dr) , (B.l) 

where / and g satisfy the integrability conditions for u and v, 

d 2 u _ d*v_ (Pv_ _ d 2 v 
dtdr drdt 1 dtdr drdt 

Without loss of generality, we shall assume that they are all strictly positive, 

/ > 0, g > 0. (B.3) 

From the above we can see that the rays moving along the hypersurfaces u = Const, are outgoing, while the ones 
moving along the hypersurfaces v = Const, are ingoing. Then, it is easy to show that, in terms of u and v, the metric 
(A.l) takes the form 

ds 2 = I 2 ^2e 2l7 ^dudv - R 2 (u,v)d8 2 } , (B.4) 
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where 



a(u,v) = -hn(2fg), R(u,v)=rS. (B.5) 



On the other hand, from Eq.(B.l) we find that 



df_ _ J_ _$ &t_ _ 1_ _ $ 
9m ~ 2/ 6 ' 0« ~ 2<? 6 ' 
cV 1 _ 9 dr 1 _^ 

^ = -27 e ' ar^ 6 • (R6) 



Then, the expansions of the outgoing and ingoing null geodesies are defined as [11], 



PR PR 

6 n = V x n x = e- 2 -|| = ^ (e-*fl, t - e~* R, r ) , (B.7) 



where Va denotes the covariant derivative, and l x (n x ) is the null vector defined the outgoing (ingoing) null geodesies 
and given by 

l *=d& = *> n * s &* = * (R8) 

For the self-similar solutions of the zeroth kind, Eq.(B.7) takes the form, 

Oi = f r {(i + y)e-*-ye^-*}, 

e n = -jt- r {{l + y)e-* ■+ye^-"}, (a = 0), (B.9) 
while for the ones of the first or second kind, it becomes 

°' = 4r \ a (! + V) e ~* + ye e+ ^- 1 ^ a -*\ , 

9 n = -JL-.{a{l + y)e- 9 - y e°+(°-W"-*} , (a ? 0). (B.10) 
Setting ; = in the above expressions, we shall find the location of the apparent horizons. 
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